Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory 
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Neutron matter presents a unique system for chiral effective field theory (EFT), because all 
many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N 3 LO). 
We present the first complete N 3 LO calculation of the neutron matter energy. This includes the 
subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. 
We find relatively large contributions from N 3 LO 3N forces. Our results provide constraints for 
neutron-rich matter in astrophysics with controlled theoretical uncertainties. 
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The physics of neutron matter ranges from universal 
properties at low densities to the structure of extreme 
neutron-rich nuclei and the densest matter we know to 
exist in neutron stars. For these extreme conditions, con- 
trolled calculations with theoretical error estimates are 
essential. Chiral EFT provides such a systematic expan- 
sion for nuclear forces [l[. This is particularly exciting 
for neutron matter and neutron-rich systems, because all 
three- and four-neutron forces are predicted to N 3 LO Q. 

Neutron matter based on chiral EFT has been studied 
using lattice simulations Q at low densities, n < no/10 
(with saturation density no = 0.16 fm~ 3 ), and following 
an in- medium chiral perturbation theory approach 0, [B| , 
where low-energy couplings are adjusted to empirical nu- 
clear matter properties. In addition, the renormalization 
group (RG) has been used to evolve chiral EFT interac- 
tions to low momenta Q , which has enabled perturbative 
calculations for nucleonic matter 0, 0] ■ While these con- 
strain the properties of neutron-rich matter to a much 
higher degree than is reflected in neutron star model- 
ing 0, the dominant uncertainties are due to 3N forces, 
which were included only to N 2 LO. A consistent inclusion 
of higher-order many-body forces is therefore key. 

Here we present the first calculations at nuclear den- 
sities based directly on chiral EFT interactions without 
RG evolution. To this end, we have studied the per- 
turbative convergence of chiral two-nucleon (NN) poten- 
tials for neutron matter in detail, and found that the 
available N 2 LO and N 3 LO potentials with lower cutoffs 
A = 450 — 500 MeV are perturbative. This is supported 
by small Weinberg eigenvalues at low energies indicat- 
ing the perturbative convergence in the particle-particle 
channel [6]. In neutron matter, it comes as a result of 
effective range effects which weaken NN interactions 
at higher momenta, combined with weaker tensor forces 
among neutrons, and with limited phase space at finite 
density due to Pauli blocking fiol j . 

At the NN level we use the N 2 LO and N 3 LO po- 
tentials developed by Epelbaum, Glockle and Meifiner 
(EGM) [IH with A/A = 450/500 and 450/700 MeV (A/ A 
denotes the cutoff in the Lippmann-Schwinger equation 



and in the two-pion-exchange spectral-function regular- 
ization, respectively). We also use the A = 500 MeV 



N 3 LO NN potential of Entem and Machleidt (EM) 12], 



which is most commonly used in nuclear structure calcu- 
lations. The larger A = 550 - 600 MeV NN potentials of 
EGM and EM have been found to be nonperturbative [HI 
and are therefore not included. Moreover, the LO NN 
contact couplings in the 600/600 and 600/700 EGM po- 
tentials break Wigner symmetry perturbatively (at the 
interaction level), with a repulsive spin- independent Cs 
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FIG. 1. (Color online) Neutron matter energy per particle 
as a function of density including NN, 3N and 4N forces at 
N 3 LO. The three overlapping bands are labeled by the differ- 
ent NN potentials and include uncertainty estimates due to 
the many-body calculation, the low-energy a constants and 
by varying the 3N/4N cutoffs (see text for details). For com- 
parison, results are shown at low densities (see also the inset) 
from NLO lattice 0] and Quantum Monte Carlo (QMC) sim- 
ulations [2^] , and at nuclear densities from variational (APR; 
the different points are with/without boost corrections) 
and Auxiliary Field Diffusion MC calculations (GCR) 
based on adjusted nuclear force models. 
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and an unnaturally large spin-dependent Ct ~ Cg, lead- 
ing to unexpectedly large CV-dependent 3N forces. 

In this Letter, we include for the first time all N 3 LO 3N 
and 4N forces, which have been derived only recently 0- 
17| . in addition to the N 2 LO 3N forces. Figure [T] shows 
our complete N 3 LO calculation of the neutron matter 
energy as our main result, where the bands include esti- 
mates of the theoretical uncertainties due to the many- 
body calculation and in the many-body forces. 

For neutrons, only the two-pion-exchange 3N forces 
contribute at N 2 LO For the corresponding low- 

energy constants c± and C3, we take the range of val- 
ues from a high-order analysis (lij . at N 2 LO: c\ = 
-(0.37 - 0.81) GeV- 1 and c 3 = -(2.71 - 3.40) GeV -1 
(which includes the a values in the EGM and EM NN 
potentials), and when the N 2 LO 3N forces are included 
in an N 3 LO calculation: c x = -(0.75 - 1.13) GeV" 1 and 
c 3 = -(4.77- 5.51) GeV -1 . It has been shown @ that 
the N 2 LO 3N force contributions in neutron matter can 
be to a good approximation calculated at the Hartree- 
Fock level. In this first calculation, we therefore evaluate 
the N 3 LO 3N and 4N force contributions to the energy 
per particle E/N at the Hartree-Fock level. The A-hody 
contributions arc then given by 
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with short-hand notation i = V.i<7i. A a denotes the A- 
body antisymmetrizer and = #(/cf — hi) the Fermi- 
Dirac distributions at zero temperature. We use a 
Jacobi-momenta regulator; in terms of ki given by fa = 
exp[-((fc 2 + . . .+fc 2 -ki • k 2 -. . .-k A _!-k A )/(AA 2 ])«-p] 
with n cxp = 4 and 3N/4N cutoff A = 2 - 2.5 fm \ For 
the nucleon and pion mass, we use m — 938.92 MeV and 
rn^ = 138.04 MeV, and for the axial coupling g A = 1.29 
and the pion decay constant /„. — 92.4 MeV. 
Chiral 3N forces at N 3 LO can be grouped into 
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where we take the long-range parts, the subleading 
two-pion-exchange, the two-pion-one-pion-exchange and 
the pion-ring 3N forces, from Ref. [15j, and the short- 
range parts, the two-pion-exchange-contact and relativis- 
tic 1/m-corrections 3N forces from Ref. 16]. In Fig. [5J 
we give the individual Hartree-Fock contributions to the 
neutron matter energy. The evaluation is aided because 
parts of the different 3N force topologies vanish for neu- 
trons, and the results have been checked by two indepen- 
dent calculations. The details of the calculation will be 
presented in a future paper. At the Hartree-Fock level, 
the 3N/4N contributions change by < 5% if the cutoff is 



taken to infinity (i.e., Jr = 1), but we will also include 
N 2 LO 3N forces beyond Hartree-Fock. This requires a 
consistently used regulator. Estimates of the theoretical 
uncertainty are provided by varying the 3N/4N cutoff. 

The two-pion-exchange 3N forces at N 3 LO can be 
largely written as shifts of the low-energy constants, 



8ci = -0.13 GeV -1 and 5c 3 = 0.89GeV _i [15J] of the 
N 2 LO 3N forces, plus a smaller contribution. The result- 
ing energy of about —1.5 MeV per particle at saturation 
density no in Fig. [5] is ~ 1/3 of the N 2 LO 3N energy, as 
expected based on the chiral EFT power counting. In 
contrast, the two-pion-one-pion-exchange 3N force con- 
tributions, which include 14 diagrams, are relatively large 
with —3.6 MeV per particle at saturation density. Of sim- 
ilar, but opposite size are the pion-ring 3N force contri- 
butions, with +3.3 MeV per particle at no- The shorter- 
range parts of N 3 LO 3N forces depend on the momentum- 
independent NN contacts, Ct and Cg, which we take 
consistently from the N 3 LO EM/EGM potential used. 
The contributions from the two-pion-exchange-contact 
3N forces include 11 diagrams and depend only on Ct- 
The resulting energy ranges from —2.8 to +1.3 MeV at no 
depending on the NN potential used. These larger 3N re- 
sults at N 3 LO are consistent with contributions from the 
large & constants at N 4 LO exactly in these three topolo- 
gies 18j . This shows that higher-order many-body forces 



still need to be investigated and that a chiral EFT with 
explicit A excitations may be more efficient, since this 
would capture these effects already at N 3 LO. Finally, the 
relativistic-corrections 3N forces depend also on fi$ and 
fig [HI and contribute at the few hundred keV level. 

The 4N force contributions in Fig. [2] are an order of 
magnitude smaller than those from the N 3 LO 3N forces 
and of similar size as the 3N relativistic corrections. We 
follow the 4N force notation V a through V n of Ref. [13], 
and include the direct and all 23 exchange terms. Due 
to the spin-isospin structure, only 3 topologies contribute 
to neutron matter: the three-pion-exchange 4N forces V a 
and V e , and the pion-pion- interaction 4N forces . The 
4N forces V k and V n involving the contact Ct vanish in 
neutron matter due to their spin structure. We find a to- 
tal 4N force contribution of —174+ 10 keV per particle at 
no- The V e and V? energies largely cancel [l9j, and their 
sum agrees with the very small ~ — 20keV per particle 
at no of Ref. [13] , which considered these two parts. 

Since diagrams beyond Hartree-Fock involving NN in- 
teractions and N 2 LO 3N forces (in particular with the 
larger Ci at N 3 LO 3N and without RG evolution) pro- 
vide non-negligible contributions we include all such 
diagrams to second order, as well as particle-particle dia- 
grams to third order, which is technically possible based 
on Ref. ■ In addition to using NN potentials with dif- 
ferent cutoffs and varying the 3N/4N cutoffs, we include 
estimates of the theoretical uncertainties of the Cj con- 
stants and in the convergence of the many-body calcu- 
lation. The latter is probed by studying the sensitiv- 
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FIG. 2. (Color online) Energy per particle versus density for all individual N 3 LO 3N and 4N force contributions to neutron 
matter at the Hartree-Fock level. The bands are obtained by varying the 3N/4N cutoff A = 2 — 2.5 fm - . For the two-pion- 
exchange-contact and the relativistic-corrections 3N forces, the different bands correspond to the different NN contacts, Ct 
and Cs, determined consistently for the N 3 LO EM/EGM potentials. The inset diagram illustrates the 3N/4N force topology. 



ity of the energy to the single-particle spectrum used. 
We find that the energy changes from second to third 
order, employing a free or Hartree-Fock spectrum, by 
0.8, 0.4, 1.3 MeV (1.4, 0.9, 2.7 MeV) per particle at n /2 
(n ) for the EGM 450/500, 450/700, EM 500 N 3 LO po- 
tentials, respectively. The results, which include all these 
uncertainties, are displayed by the bands in Fig. [TJ Un- 
derstanding the cutoff dependence and developing im- 
proved power counting schemes remain important open 
problems in chiral EFT [2l[ . For the neutron matter en- 
ergy at no, our first complete N 3 LO calculation yields 
14.1 — 21.0 MeV per particle. If we were to omit the 
results based on the EM 500 N 3 LO potential, as it con- 
verges slowest at no, the range would be 14.1 — 18.4 MeV. 

As we find relatively large contributions from N 3 LO 
3N forces, it is important to study the EFT convergence 
from N 2 LO to N 3 LO. This is shown in Fig. [3] for the 
EGM potentials (N 2 LO is not available for EM), where 
the N 3 LO results are found to overlap with the N 2 LO 
band across a ±1.5 MeV range around 17 MeV at satura- 
tion density. As expected from the net-attractive N 3 LO 
3N contributions in Fig. [2j the N 3 LO band yields lower 
energies. For the N 2 LO band, we have estimated the the- 
oretical uncertainties in the same way, and the neutron 
matter energy ranges from 15.5 — 21.4 MeV per particle 



at no. The theoretical uncertainty is reduced from N 2 LO 
to N 3 LO to 14.1 - 18.4 MeV, but not by a factor ~ 1/3 
based on the power counting estimate. This reflects the 
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FIG. 3. (Color online) Neutron matter energy per particle as 
a function of density at N 2 LO (upper/blue band that extends 
to the dashed line) and N 3 LO (lower/red band). The bands 
are based on the EGM NN potentials and include uncertainty 
estimates as in Fig. [T] 



4 



> 



20 



N 3 LO 
(this work) 




Tamii et al. (2011) 



Hebeler efal. (2010) 



30 



32 

S v [MeV] 



34 



36 



FIG. 4. (Color online) Range for the symmetry energy S v 
and its density dependence L obtained at N 3 LO (this work) 
versus including 3N forces at N 2 LO (Hebeler et al. [8]). For 
comparison 25}], we show constraints obtained from energy- 
density functionals for nuclear masses (Kortelainen et al. [2611') 
and from the 208 Pb dipole polarizability (Tamii et al. [271]). 

large Cj 3N contributions at N 4 LO, and is similar to the 
convergence pattern observed in chiral NN potentials [l| . 

The neutron matter energy in Fig. [T] is in very good 
agreement with NLO lattice results [3J and Quantum 
Monte Carlo simulations [22| at very low densities (see 
also the inset) and approximately reproduces the scaling 

~ 0.5 ig^j, which we attribute to effective-range effects 
combined with low cutoffs [9|. At nuclear densities, we 
compare our N 3 LO results with variational calculations 
based on phenomenological potentials (APR) |23[ , which 
are within the N 3 LO band, but do not provide theoret- 
ical uncertainties. In addition, we compare the density 
dependence with results from Auxiliary Field Diffusion 
MC calculations (GCR) [13] based on nuclear force mod- 
els adjusted to an energy difference of 32 MeV between 
neutron matter and the empirical saturation point. The 
density dependence is similar to the N 3 LO band, but the 
GCR results are higher below 0.05 fnV 3 . 

The N 3 LO band provides key constraints for the nu- 
clear equation of state and for astrophysics. Figure 2] 
shows, following Ref. 25[, the allowed range for the 
symmetry energy S v and its density dependence L — 
3nod n S v (no) (for details on the determination of S v 
and L see Ref. 0). Compared to the results from 
RG-evolved chiral interactions with 3N forces at N 2 LO 
only [8j, we find the same correlation (with the same 
slope), but not as tight due to the additional density de- 
pendences at N 3 LO. The N 3 LO ranges for S v and L are 
S v = 28.9 - 34.9 MeV and L = 43.0 - 66.6 MeV. The 
two neutron-matter bands in Fig. [?] are complementary, 
because the RG evolution in Hebeler et al. [8j improves 



the many-body convergence, while the band presented in 
this work is the first consistent N 3 LO calculation. The 
predicted N 3 LO range, as well as that of Hebeler et al. [H, 
are in agreement with constraints obtained from energy- 
density functionals for nuclear masses [26| and from the 
208 Pb dipole polarizability j^. In the future, the N 3 LO 
band can be narrowed further by a higher-order many- 
body calculation with N 3 LO 3N forces and by taking into 
account A excitations (explicitly or through large Cj con- 
tributions at N 4 LO Il8[). Combined with the heaviest 
2Mq neutron star [28| and a general extension to high 
densities 0, our N 3 LO energy band leads to a radius 
range of R = 9.7 — 13.9 km for a typical 1.4M neu- 
tron star, in remarkable agreement with Ref. f8|. For an 
alternative determination using in-medium chiral pertur- 
bation theory for all densities see Ref. [f|. 

We have presented the first complete N 3 LO calculation 
of the neutron matter energy, including NN, 3N and 4N 
forces, with the first application of N 3 LO 3N forces to 
many-body systems. The significant contributions from 
N 3 LO 3N forces show that their inclusion will also be 
very important for nuclear structure and reactions. Our 
results provide constraints for the nuclear equation of 
state and for neutron-rich matter in astrophysics, and 
highlight the exciting role neutron matter and neutron- 
rich systems play in chiral EFT, where all many-neutron 
forces are predicted. The large contributions from N 3 LO 
3N forces signal the importance of A contributions at 
nuclear densities. 
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